Geometric methods for the most general Ginzburg-Landau model with two order 

parameters 
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The Landau potential in the general Ginzburg-Landau theory with two order parameters and all 
possible quadratic and quartic terms cannot be minimized with the straightforward algebra. Here, 
a geometric approach is presented that circumvents this computational difficulty and allows one to 
get insight into many properties of the model in the mean-field approximation. 



Introduction. The Ginzburg-Landau (GL) theory of- 
fers a remarkably economic description of phase tran- 
sitions associated with breaking of some symmetry 
This breaking is described with an order parameter ip: 
the high symmetry phase corresponds to ^ = Q, while 
the low symmetry phase is described by ^/^ 7^ 0. In order 
to find when a given system is in high or low symmetry 
phase, one constructs a Landau potential that depends 
on the order parameter, and then find its minimum. Its 
classic form is 



(1) 



Higher order terms odV'l'*) are usually assumed to be neg- 
ligible. The values of the coefficients a and h and their 
dependence on temperature, pressure, etc. can be either 
calculated from a microscopic theory, if it is available, 
or considered as free parameters in a phenomenological 
approach. The phase transition associated with the sym- 
metry breaking takes place when an initially negative a 
becomes positive, and the minimum of the potential ([T]) 
shifts from zero to 1-01 — a/b. 

Many systems are known, in which two competing or- 
der parameters (OP) coexist. Among them are general 
0(m) © 0(n)-symmetric models, |2|, models with two in- 
teracting A^- vector OPs with 0{N) symmetry, [3]; spin- 
density- waves in cuprates, 5*0(5) theory of antiferro- 
magnetism and superconductivity, [5|; multicomponent, 
[ill, spin-triplet p- wave, Q, and two-gap, supercon- 
ductivity, with its application to magnetism in neutron 
stars, [10|; two-band superfluidity, [IjJ; various mech- 
anisms of spontaneous breaking of the electroweak sym- 
metry beyond the Standard Model such as the two-Higgs- 
doublet model (2HDM), jHI. 

To describe such a situation within the GL model, 
one constructs a Landau potential similar to ([T]), which 
depends on two order parameters, ^1 and 'ip2- Coeffi- 
cients of this potential, and bi , can be considered in- 
dependent, although in each particular application they 
might obey specific relations. One thus arrives at the 
most general two-order-parameter (20P) GL model with 
quadratic and quartic terms. 

A natural question arises: what is the ground state 
of the most general 20P GL model? A rather surpris- 



ing fact is that this question cannot be answered by a 
straightforward calculation. Differentiating the Landau 
potential in respect to ipi leads to a system of coupled 
algebraic equations that cannot be solved explicitly. 

In this Letter we show that despite this computational 
difhculty, one can still learn much about the most general 
20P GL model. Namely, one can study the number and 
the properties of the minima of the Landau potential, 
classify possible symmetries and study when and how 
they are broken. In short, one can describe the phase 
diagram of the model, at least in the mean-field approx- 
imation, without explicitly minimizing the potential. 

There exists, in fact, an extensive literature dating 
back to 1970's on minimization of group-invariant po- 
tentials with several OPs with the aid of stratification of 
the orbit space, see e.g. [H Here we show that in 
the case of two order parameters realizing the same group 
representation the analysis can be extended much farther 
than in the general case, with important physical conse- 
quences. For a particular application of this formalism 
to the 2HDM, see 
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The formalism. Let us focus on the simplest case when 
two OPs ipi and ip2 are just complex numbers. The most 
general quadratic plus quartic Landau potential is 



r(^r^2) + 651^1 



■ 63|'0lHV'2 



(2) 



c.c. 



It contains 13 free parameters: real oi, 02, 61, 62, ^3 and 
complex as, 64, 65, bg. For the illustration of the main 
idea, we place no restriction on jV'il from above. Note 
that potential ([2]) contains quartic terms that mix tpi and 
"02 • Such terms are usually absent in particular applica- 
tions of the 20P GL model (for a rare exception, see [i|), 
but in the approach presented here it is essential that all 
possible terms are included from the very beginning. 

Once potential ^ is written, the physical nature of 
OPs becomes irrelevant. One can consider them as com- 
ponents of a single complex 2-vector $ = (V'i,V'2)"^- 
The key observation is that the most general poten- 
tial jS]) keeps its generic form under any regular linear 
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transformation between tpi and tp2' $ ^ = T • $, 
T € GL{2,C). It can be also accompanied with a suit- 
able transformation of the coefficients a^, 6^, so that one 
arrives at exactly the same potential as before. Thus, 
the problem has some reparametrization freedom with 
the reparametrization group GL{2,C). Among 13 free 
parameters, only 6 play crucial role in shaping the phase 
diagram of the model, while the other 7 just reflect the 
way we look at it. 

Let us now introduce a four- vector r'^ — (ro, r^) — 
(cf>V^$) with components 



ro = 




2Re 

2Im iri^2) 



(3) 



IV'2 



Here, index fx = 0,1,2,3 refers to components in the 
internal space and has no relation with the space-time. 
Multiplying ipi by a common phase factor does not 
change r^, so each parametrizes a U (l)-orbit in the ipi- 
space. Since the Landau potential is also ?7(l)-invariant, 
it can be defined in this 1 -f- 3-dimensional orbit space. 

The SL{2, C) C GL{2, C) group of transformations of 
$ induces the proper Lorentz group 5*0(1, 3) of transfor- 
mations of r^. This group includes 3D rotations of the 
vector Vi as well as "boosts" that mix tq and r^, so the 
orbit space gets naturally equipped with the Minkowski 
space structure. Since rg > and r^r^ = ~ J2 = 0, 
the orbit space is given by the "forward lightcone" LC+ 
in the Minkowski space. 

All this allows us to rewrite ([2]) in a very compact form: 



(4) 



with 



= - (ai -I- a2, -2Re as, 2Im as, -oi -I- 02) , 



Finding eigenvalues B'^'^ explicitly in terms of a^, bi 
requires solution of a fourth-order characteristic equa- 
tion, which constitutes one of the computational difhcul- 
ties of the straightforward algebra. We reiterate that in 
our analysis we never use these explicit expressions. The 
analysis relies only on the fact that the eigenvalues are 
real and satisfy ([6]). 

Minima of the Landau potential. Let us first find how 
many extrema the potential Q can have in the orbit 
space. Since extrema lie on the surface of LC^ , we use 
the Lagrange multiplier method to arrive at the following 
simultaneous equations: 



B^,{r^)-\{T^)=A^, (r^)(r^) = 0. 



(7) 



Here, (r^) labels the position of an extremum. This sys- 
tem cannot be solved explicitly in the most general case, 
however one can establish how many extrema a given po- 
tential has. 

To find it, we rewrite (r^) — {ro)ni, where Ui is a unit 
3D vector; then ^ becomes 



[Aa-{Ba-B,){ra)]n,^A^ 



(8) 



Assume for simplicity that Aq > and Bi < B2 < 
B3. The l.h.s. of ^ at fixed (rg) and all unit vec- 
tors Hi parametrizes an ellipsoid with semiaxes Aq — 
{Bo — Bi){ro). As (ro) increases from zero to infin- 
ity, this ellipsoid first shrinks, then grows, collapsing at 
Tq*' = Aq/(Bq — Bi) to planar ellipses. One can see that 
during these transformations for (rp) ~ [0,00) it sweeps 
at least once the entire Minkowski space and at least 
twice the interior of the sphere of radius Aq. In ad- 
dition, there are two cusped regions, such as shown in 
Fig. [U whose interior is swept twice more. So, by check- 
ing whether Ai lies inside these regions, one can get the 
number of solutions of ^ without finding them explic- 
itly. 
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Here, b^^ = 61 ± 62, 655 = 65 ± fog- 

One usually requires that the quartic term of the po- 
tential increases in all directions in the OP-space. In the 
orbit space, this was proved in [14] to be equivalent to 
the statement that is diagonalizable by an 50(1,3) 
transformation and after diagonalization it takes form 
Bf,i, = diag(Bo, -Bi, -B2, -S3) with 



Bo > Bi,B2, i?s 



(6) 



Since ' 



= 0, the matrices B^"' and B^"' = B^"'~Cg^"' 



are equivalent. This degree of freedom in the definition 
of B^^'^ shifts all the eigenvalues by a common constant. 




FIG. 1: The envelope of ellipsoids for Tq^^ < ro < 1 q 
j4i-space. 



rn^^ in the 



In the 1 -I- 3-dimcnsional space of yl^ , these 3D regions 
serve as bases of corresponding conical regions with dif- 
ferent numbers of extrema of the potential. Namely, at 
least one non-trivial solution of ([7]) exists, if A^ lies out- 
side the past lightcone LC~ (otherwise the global mini- 
mum of the potential is at the origin, {ipi) =0). If lies 
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inside the future lightcone LC"*", then there are at least 
two non- trivial solutions. If in addition lies inside 
one or both caustic cones, then two additional extrema 
per cone exist. In total, the potential can have up to six 
non-trivial extrema in the orbit space. This fact was also 
found independently in (l6j . 

The above construction does not distinguish a mini- 
mum from a saddle point (with condition ([6]), there are no 
non-trivial maxima in the orbit space) . A straightforward 
method for finding a minimum, which consists in check- 
ing that the hessian eigenvalues are all non-negative, is 
again of little use here. Instead, one can still use geome- 
try to study the properties of the minima. 

As described above, physically realizable points of the 
orbit space lie on LC'^. Nevertheless, let us consider ex- 
pression in the entire Minkowski space. Let us define 
an equipotential surface as a set of all vectors r'' 
with the same value of V. These equipotential surfaces 
do not intersect, are nested into each other, and have very 
simple geometry: they are second-order 3-dimensional 
manifolds (3-quadrics). Since LC^ is also a specific 3- 
quadric, finding points in the orbit space with the same 
value of V amounts to finding intersections of these two 
3-quadrics. In particular, to find a local minimum of the 
potential in the orbit space, one has to find an equipoten- 
tial surface that touches LC'^ in an isolated point (we say 
that two surfaces "touch" if they have parallel normals 
at the intersection points). The global minimum cor- 
responds to the unique equipotential surface that only 
touches but never intersects iC+. The fact that the 
search for the global minimum is reformulated in terms 
of contact of two 3-quadrics leads to several Propositions 
listed below. 

Let us now find how many among the extrema are 
minima. Let us fix B'^'^ and move A'^ continuously in 
the parameter space. We first note that the signature of 
the hessian can change only when the total number of 
extrema changes. A saddle point cannot simply become 
a minimum; it can only split into several extrema, one 
of them being a minimum, or it can merge with other 
extrema to produce a minimum. Therefore, the conical 
3-surfaces described above {LC^ , LC^, caustic cones) 
separate regions in the A'^-space with a definite number 
of minima. One can then check a representative point 
— {Aq, 0, 0, 0) (in the basis where i?^'' is diagonal) 
of the innermost region in the A'^ space and find that 
there are two distinct minima in this case. This proves 
the following Proposition: 

Proposition 1. The most general quadratic plus quartic 
potential with two order parameters can have at most two 
distinct local minima in the orbit space. 

Symmetries and their violation. The potential can 
have an additional explicit symmetry, i.e. it can remain 
invariant under some transformations of ipi (or coeffi- 
cients) alone. If the position of the global minimum is 



also invariant under the same group of transformations, 
we say that the symmetry is preserved; otherwise, we 
say that the explicit symmetry is spontaneously violated. 
In most applications, the Landau potential does possess 
some explicit symmetry, so whether it is preserved or vi- 
olated can have profound physical consequences. 

An explicit symmetry corresponds to such a map of 
the Minkowski space that leaves invariant, separately, 
B^i^r'^r'^ and Af^r'^. In a general situation, it might be far 
from evident that the potential has any explicit symme- 
try. The following criterion helps recognize the presence 
of a hidden explicit symmetry and tells what symmetry 
it is: 

Proposition 2. Suppose that the potential is explicitly 
invariant under some transformations of . Let G he 
the maximal group of such transformations. Then: 
(a) G is non-trivial if and only if there exists an eigen- 
vector of B^i, orthogonal to A^j,; 

(h) G is one of the following groups: Z2, (^2)^7 , 
0(2), 0(2) X Z2, 0(3). 

The proof will be given in [l3| (see also [l^ for a similar 
statement in 2HDM). In the case of a discrete symmetry 
the following Propositions can be easily proved: 

Proposition 3. The maximal violation of any discrete 
explicit symmetry consists in removing only one Zi fac- 
tor: {Z2f^{Z2f-\ 

Proposition 4. For any explicit discrete symmetry, 
minima that preserve and spontaneously violate this sym- 
metry cannot coexist. 

Both Propositions follow from Proposition [1] and the 
fact that the set of all minima preserves the explicit sym- 
metry group G. 

If a discrete symmetry is spontaneously violated, then 
there are two generate minima in the orbit space. One 
can also prove the converse, i.e. the two degenerate min- 
ima can arise only via spontaneous violation of a discrete 
symmetry of the potential, [l5| . Thus, the criterion for 
the spontaneous violation of a discrete symmetry is that 
A^ lies inside a caustic cone associated with the largest 
eigenvalue. 

For a concrete example, suppose that all eigenvalues 
of B^^ are distinct and that A3 = 0, while other compo- 
nents Ai 7^ 0. This potential has an explicit Z2 symme- 
try generated by reflections of the third coordinate. The 
global minimum spontaneously violates this symmetry 
(i.e. (7-3)7^0), if B3>Bi,B2 and 



Al 



< 



^0 



{B3 - Bi)2 {B3 - B2Y {Bo - B3) 



(9) 



The proof is based on the "shrinking ellipsoid" construc- 
tion described above and will be given in detail in [l3| . 

Local order parameters. In this Letter we have illus- 
trated the idea using the global OPs ipi- The approach 
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can be easily extended to models, where OPs are defined 
locally, ipi{x)- In this case, one considers the free-energy 
functional Fltpi] — J (Px[K{ipi) + V{ipi)] with kinetic 
term K and potential V. In the general model the ki- 
netic term must include all quadratic combinations of 
the gradient terms: 

K - Ki|^^i|2 + K2|^V2|'+^£30V'l)*(^V'2)+C.C. , (10) 

where D is either V or the covariant derivative. It can 
be rewritten in the reparametrization invariant form K = 
Kfj^pf^ with = (i3$)t(T^(z3$) and defined in terms 
of Ki similarly to defined in terms of Oj. 

The presence of Kfj_ leads only to minor complications 
of the above analysis. All the conclusions about the num- 
ber of extrema and minima remain unchanged, however 
one should now distinguish symmetries of the potential 
and of the entire free energy functional, see [17[ for de- 
tails. 

Two local OPs also lead to the existence of quasitopo- 
logical solitons. This possibiUty has been known for 
some time; for example, in 19], a soliton in the onc- 
dimensional two-band superconductor with a simple in- 
terband interaction term was described. Such a soliton 
corresponds to the relative phase between the two con- 
densates that changes continuously from zero to 27r as 
X goes from — oo to +oo, and it is stable against small 
perturbations. The general origin of such quasitopologi- 
cal solitons is obvious from the above construction. The 
orbit space of all non-zero configurations of OPs is the 
forward lightcone without the apex, which is homotopi- 
cally equivalent to a 2-sphcre . Depending on the exact 
shape of the potential, it can support either closed lin- 
ear paths, which correspond to domain walls, or closed 
2-manifolds, which corresponds to strings. 

MuUicomponent order parameters. In many physical 
situations one encounters multicomponent OPs. Exam- 
ples include 2IIDM, superfluidity in ■^He, spin-density 
waves, etc. The formalism developed here works also 
for these cases. Just to mention some characteris- 
tic features, leaving a detailed discussion for [l3l, we 
note that 5?7(A^)-symmetric potential depends on A^- 
vectors only via combinations (ipjipj)- Since in general 

(■01 '02) (^2 ^i) 7^ |'0iHV'2p, one gets a new term in the 
potential ^ with its own coefficient. The definition of 
r^ remains the same, but r'^r^ > 0, so r'^ can lie not 
only on the surface, but also in the interior of LC~^. This 
makes definition of i?^^ unique, fixes its eigenvalues, and 
depending on their signs, one has to consider separately 
several cases. Modifications to the overall results are mi- 
nor, see IJ, |l5| for a 2HDM analysis. One obtains a 
new phase, with (r'^) lying inside LC^ , which has differ- 



ent symmetry properties (in 2HDM it corresponds to a 
charge-breaking vacuum). One can easily formulate con- 
ditions when it is the global minimum of the model, so 
the phase diagram remains equally tractable in this case. 

In conclusion, we considered the most general 
Ginzburg-Landau model with two order parameters, in- 
cluding all possible quadratic and quartic terms in the 
Landau potential. Since the minimization of the poten- 
tial cannot be done explicitly, we developed a geomet- 
ric approach based on the reparametrization properties 
of the model and used it to study the ground state of 
the model. Future research should include dynamics of 
the fluctuations above the ground state, corrections to 
the potential beyond the quartic term, renormalization 
group flow, as well as modiflcations of the results at fl- 
nite temperature and in the presence of external fields. 
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